Direct Measurement of the Out-of-Plane Spin Texture in the Dirac Cone Surface 

State of a Topological Insulator 
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We have performed spin- and angle-resolved photoemission spectroscopy of Bi2Te3 and present 
the first direct evidence for the existence of the out-of-plane spin component on the surface state 
of a topological insulator. We found that the magnitude of the out-of-plane spin polarization on a 
hexagonally deformed Fermi surface of Bi2Te3 reaches maximally 25% of the in-plane counterpart 
while such a sizable out-of-plane spin component does not exist in the more circular Fermi surface 
of TlBiSe2, indicating that the hexagonal deformation of the Fermi surface is responsible for the 
deviation from the ideal helical spin texture. The observed out-of-plane polarization is much smaller 
than that expected from existing theory, suggesting that an additional ingredient is necessary for 
correctly understanding the surface spin polarization in Bi2Te3. 

PACS numbers: 73.20.-r, 79.60.-i, 71.20.-b, 75.70.Tj 



The topological insulators (TIs) materialize a new 
state of quantum matter where an unusual gapless metal- 
lic state appears at the edge or the surface of a band in- 
sulator due to a topological principle. The surface state 
(SS) of three-dimensional TIs is characterized by a Dirac- 
cone dispersion which has been shown to have a helical 
spin structure where the spin vector points parallel to 
the surface and perpendicular to the momentum k, as 
shown in Fig. 1(a) IH [I]- Because of this helicity in the 
spin direction and the protection by the time-reversal 
symmetry, the Dirac fermions in the TIs are immune to 
the backward scattering |3| and are not very sensi- 
tive to nonmagnetic impurities or disorder. This peculiar 
situation provides a platform for novel topological phe- 
nomena such as the emergence of Majorana fermions in 
the proximity-induced superconducting state, and indeed 
many theoretical models or experimental results were de- 
veloped or interpreted relying essentially on this simple 
helical spin structure [H IH, IJ However, it is unclear at 
the moment to what extent such a simple spin texture 
is adequate. In fact, recent theoretical studies predicted 
that when the Fermi surface (FS) of the surface Dirac 
state is hexagonally deformed, the spin structure starts 
to obtain a finite out-of-plane (OP) component It 
was also predicted 0, Q that such a component can be 
as large as the in-plane (IP) counterpart in Bi2Te3 which 
shows the strongest hexagonal FS warping among known 
TIs 0,11. The deviation from the simple helical spin tex- 
ture would be a key to better understanding the peculiar 
SS of the TIs, whereas an experimental verification has 
not yet been successfully made. 

In this Letter, we demonstrate that such an OP spin 
component is indeed present in Bi2Te3, by determining 
for the first time the detailed IP and OP spin texture 
at various k locations of the Dirac-cone SS using the 



low-energy, spin- and aiigle-resolved photoemission spec- 
troscopy (SR-ARPES) [lai. We also elucidated that the 
magnitude of the OP spin polarization is related to the 
strength of the hexagonal warping of the FS. We discuss 
the present result in relation to theoretical proposals as 
well as its implications for potential applications. 

High-quality single-crystals of Bi2Te3 were grown by 
melting nearly stoichiometric mixtures of high-purity 
(99.9999% pure) elemental shots of Bi and Tc in a quartz 
glass tube at 850°C for 48 h and slowly cooling to 500°C 
over 48 h. To tune the chemical potential to be located 
close to the bottom of the conduction band (CB), which 
was crucial for the success of this work, we used a slightly 
Bi-rich composition of Bi:Te = 0.34:0.66. ARPES mea- 
surements were performed by u sing an SR-ARPES spec- 
trometer at Tohoku University [10|, [ill . We used the Xe 
I photons (8.437 cV) to excite photoelectrons. Samples 
were cleaved in-situ along the (111) crystal plane in an 
ultrahigh vacuum of 5xl0~^^ Torr. The energy resolu- 
tions during the SR- and regular ARPES measurements 
were 40 and 6 meV, respectively. We used the Sherman 
function value of 0.07 to obtain the SR-ARPES data. 

Figure 1(b) shows the FS mapping of an n-type Bi2Te3 
sample. We clearly find a FS showing a strong hexagonal 
deformation as in previous reports (5|, |9| . This snowflake- 
like FS originates from an electronlikc surface band [Fig. 
1(c)] corresponding to the upper branch of the Dirac 
cone. The Dirac point is situated at the binding energy 
Eb of 0.25 eV, suggesting that the bottom of the bulk CB 
is located at slightly above the Fermi level (c/. Ref. P); 
in other words, the chemical potential (/i) lies slightly 
below the bottom of the bulk CB. This is consistent with 
the slightly n-type character of the crystal indicated by 
the Seebeck coefficient and is suited for examining the OP 
spin component because of the large hexagonal warping 
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FIG. 1: (Color online) (a) Schematic picture of the Dirac-cone 
SS with a simple helical spin texture, (b) ARPES intensity at 
Ef of an n-type Bi2Te3 crystal around the T point plotted as 
a function of two-dimensional wave vector k, highlighting the 
snowflakelike FS. The ARPES intensity is integrated within 
±20 meV with respect to Ef and folded by taking into account 
the crystal symmetry. Definition of the FS angle <j> and the IP 
up-spin vector (</> up) are also indicated, (c) Near-iJp ARPES 
intensity measured along the TM cut. (d),(e) SR-EDCs for 
the IP (0) and OP (z) components, respectively, measured at 
various k points indicated by open circles in (b). (f),(g) IP 
and OP spin polarizations and Pz at points A, D, and C. 



of the FS at this energy position [§] . It should be empha- 
sized here that the absence of bulk CB electrons below 
E-p makes it possible to quantitatively determine the net 
spin polarization solely from the surface band by avoiding 
the contamination from the bulk band. 

Figure 1(d) shows spin- resolved energy distribution 
curves (SR-EDCs) of Bi2Te3 for the IP spin component 
measured at various k locations around the FS indicated 
in Fig. 1(b). We define the IP spin polarization vector 
to point perpendicular to measured k [see, e.g., the thick 
(red) arrow in Fig. 1(b)] in which the "up spin" points to 
the clockwise direction. A peak near the Fermi level (Ep ) 
in the up-spin EDCs originates from the upper branch of 
the Dirac cone, while a shoulder feature at Eb > 0.2 eV 
is a tail of the bulk valence band. In contrast to the small 
differences between the up- and down-spin EDCs for the 
shoulder feature at E^ > 0.2 eV, the difference for the 
near-i^F peak is gigantic: it is dominated by the up-spin 
component, indicating the basically spin helical nature of 
the SS where the IP spin vector points to the clockwise 
direction like in Fig. 1(a). 



The important question is whether there is a finite 
OP spin component perpendicular to the sample sur- 
face, Sz- We can directly address this issue with our SR- 
spectrometer, which allows us to simultaneously measure 
the OP component with the IP component, and Fig. 1(e) 
shows the SR-EDCs for the OP spin component. If there 
is no OP spin polarization, one should observe that the 
EDCs for Sz-np and Sz-down should completely coincide 
with each other. It turns out that this is not the case, 
and a clear Sz dependence of the EDCs was observed. 
For instance, around the Fermi vector (fcp) along the TK 
direction corresponding to the FS angle (j) of 0° (points 
A and B), the Sz-down EDC near Ep is clearly enhanced 
over to the S'z-up one around the peak top, although the 
difference between the two EDCs is much smaller than 
that in the IP polarization measurements [Fig. 1(d)]. 
Furthermore, the relative intensity of the neav-Ep peak 
for the Sz-down and S'z-up spins is reversed along an- 
other TK direction (</> = 60° , which we call here the TK' 
direction) as shown by the EDCs at point C. This demon- 
strates a sign change of the OP spin polarization within 
the SS: Spins at the k-space locations A and B point to 
the sample side, whereas those at location C point to 
the vacuum side. Notably, there is no difference between 
S'^-up and S'^-down EDCs along the TM direction (point 
D; = 30°) within our experimental accuracy, suggest- 
ing that the sign change takes place along the TM line. 
Also, one can see in Fig. 1(e) that there is no clear differ- 
ence between S'z-up and S^-down EDCs at k near the F 
point (points E and F), suggesting that the Sz component 
vanishes as the energy position on the SS approaches the 
Dirac point and the hexagonal warping weakens. 

To discuss quantitatively the present observation, we 
plot in Figs. 1(f) and (g) the energy dependence of the 
IP (P^) and OP (Pz) spin polarizations, respectively, for 
three representative k points. The spin polarization P 
is defined as P = (N^-N;)/(N^-hN;) where and N; 
are the intensity of the up- and down-spin states, respec- 
tively. One can clearly see that the P^ value increases 
as Eb approaches Ep, and it reaches ~0.6 around Ep. 
In contrast, the magnitude of |Pz| is much smaller, and 
at points A or C it reaches only ^0.15, which is about 
25% of the maximum P^ value. We have confirmed that 
this characteristic behavior is reproducible by measuring 
several samples. In addition, we have found that the ex- 
posure of an originally p-type Bi2Tc3 samples to the he- 
lium gas in the vacuum chamber leads to the emergence 
of a similar snowflakelike FS due to the surface band 
bending 12-l3| and that its spin texture is essentially 
the same as the bulk n-type samples. It is emphasized 
here that the clear observation of the OP spin compo- 
nent has become feasible by utilizing the hig h-flux and 
low-energy character of the xenon lainp 111 ] as well as 
an efficient SR- ARPES spectrometer [10[, which enable 
the high-energy or k resolution measurement with higher 
statistics. 
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To analyze the k dependence of the spin polarization 
in more detail, we plot in Figs. 2(a) and (b) the spin po- 
larizations along the x- and z-axes measured in the kx=Q 
cut [vertical solid (green) line in Fig. 1(b)] obtained by 
integrating the ARPES intensity at E^i = 0.08-0.12 eV; 
because of the rather large (40 meV) energy resolution 
of the spin-resolved measurements, the integration in this 
energy window gives a good representation of the spin po- 
larization of the surface band as a function of k. One can 
see in Fig. 2(a) that the Px value at fc^ > {ky < 0) is 
positive (negative), and its absolute value \Px\ is nearly 
constant at ~0.6 for \ky\ > 0.03 A~^. The sudden sign 
reversal in near fcy = is consistent with the IP spin 
helical texture. We also found that the z-axis spin polar- 
ization Pz shows a sign reversal along the same cut [Fig. 
2(b)], as expected from the result in Fig. 1 where the TR' 
direction at = 60° is identical to the TK' direction at 
(j> = 180° due to the three- fold rotational symmetry of 
the crystal. However, the magnitude of [P^l stays nearly 
zero at |fc.y| < 0.03 and then gradually increases 

with \ky\, reaching ~0.15 near the fcp point. This trend 
suggests that, unlike the behavior of P^;, [P^l grows sys- 
tematically as the energy position on the SS moves away 
from the Dirac point at the zone center. 

To highlight the OP spin configuration obtained from 
the present experiment, we illustrate in Fig. 2(c) the 
magnitude of Pz as a function of k. We also show numer- 
ically simulated P^ values based on the k ■ p perturbation 
theory Q which reproduces the hexagonal deformation 
of the observed FS. One can see that the measured P^ 
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FIG. 2: (color online) (a),(b) Experimental spin polarization 
along the x and z axes, Px and Pz, of the SS along the kx=Q 
cut, obtained by integrating the SR- ARPES intensity at Eb 
= 0.08-0.12 eV. (c) Schematic picture of the z component 
of the spin vector as a function of k in the SS of Bi2Te3 
determined from the present ARPES experiment. The length 
of the arrow corresponds to the magnitude of Pz ■ The data 
are symmetrized by taking into account the crystal symmetry. 
The intensity with gradual shading (red and blue) shows the 
calculated Pz value based on the theory in Ref. 0]. 



switches the sign across the f M line [straight dashed 
(green) lines] for every 60° step of (j), and its absolute 
magnitude systematically increases as one moves away 
from the T point. Such an overall experimental feature is 
in qualitative agreement with the numerical simulation. 

Next we discuss the relationship between the hexag- 
onal warping of the FS and the OP spin polarization. 
We have performed SR- ARPES measurements on the re- 
cently discovered topological insulator TlBiSe2 15-l3]. 
As shown in Fig. 3(a), the FS of TlBiSe2 is much closer 
to the circular shape as compared to the snowflakelike 
FS of Bi2Te3 [Fig. 3(b)], indicating a weaker warping ef- 
fect in TlBiSe2. While an IP spin polarization is clearly 
observed [upper traces in Fig. 3(c)], the S'z-up and Sz- 
down EDCs of TlBiSe2 overlap with each other [lower 
traces in Fig. 3(c)], indicating a negligible OP polariza- 
tion which can also be seen in the E-q dependence of Pz 
[Fig. 3(e)]. A side-by-side comparison of the Bi2Te3 data 
[Figs. 3(d) and (f)] vividly shows the difference between 
the two compounds in the magnitude of Pz- In fact, the 
ratio of the OP spin component to the IP component in 
the E-Q range of Ep-Q.2 eV is estimated to be 6±6% and 
20±4% for TlBiSe2 and Bi2Te3, respectively We have 
also confirmed the negligibly small Pz value in TlBiSe2 
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FIG. 3: (color online) (a),(b) Comparison of the FS between 
TlBiSe2 and Bi2Te3. Solid lines are the simulated FS using 
the theory in Ref. [^. (c),(d) SR-EDCs measured along the 
VK line (0 = 0°) for TlBiSe2 and Bi2Te3, respectively (e),(f) 
Corresponding Es dependences of P^. (g) Calculated Pz value 
along the kx=Q cut as a function of ky for TlBiSe2 and Bi2Te3. 
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by measuring data at various k points. These results 
point to an intimate connection between the hexagonal 
warping of the FS and the magnitude of [111 , although 
it is difficult to rule out all the other possibilities at the 
moment. To further examine the present observation, we 
reproduced the observed FS shape of TlBiSe2 and Bi2Te3 
by using the theoretical band dispersion with the hexag- 
onal warping term [see the solid curves in Figs. 3(a) 
and (b)] and calculated the value along the kx=Q cut 
[dashed lines in Fig. 3(a)] as shown in Fig. 3(g). Appar- 
ently, the calculated Pz value in TlBiSc2 is much smaller 
than that in Bi2Te3 (^0.2 and ^^0.6 at fcp, respectively), 
which is qualitatively consistent with what we observed. 

However, the magnitude of Pz is quantitatively differ- 
ent between theory and experiment, as seen in Fig. 3: 
For instance, at the k-p point along the TK direction (0 
= 0°) where \Pz\ becomes maximal on the FS, the theo- 
retical Pz of ^0.6 [Fig. 3(g)] is ~4 times larger than the 
experimental Pz of ~0.15 [Fig. 3(f)]. The smaller Pz in 
the experiment cannot be attributed to extrinsic effects 
in the experiment such as the instrumental asymmetry 
of up- and down-spin channel or a nonlinearity in the 
detector sensitivity to the actual polarization, because 
the experimental value of Pz is systematically reversed 
across ky=0 with a similar magnitude relative to ky=0 
[see Fig. 2(b)] and the experimental scattering asymme- 
try was confirmed to be proportional to the theoretical 
spin polarization by measuring Au(lll) as a reference. 
Therefore, there must be an intrinsic reason for the weak- 
ening of Pz- One possible explanation is the many-body 
effects that lead to an electronic order like the spin den- 
sity wave where the FS is nested with the wave vector 
connecting parallel segments of the hexagonally warped 
FS [3l • Another possibility is the presence of an inter- 
band scattering channel [13[ that degrades the absolute 
surface spin polarization through the scattering between 
the surface band and the non-spin-polarizcd bulk band. 
It is desirable that a concrete theoretical model is devel- 
oped to explain the unexpectedly small OP polarization. 

The presence of the OP spin polarization provides a 
novel mechanism for the "mass acquisition" of the Dirac 
particles on the surface 0], which is a prerequisite of 
exotic topological phenomena like the half-integer quan- 
tum Hall effect and the quantized magnetoelectric effect 
i^lHj. When the SS has a purely IP hchcal spin texture, 
the gapping of the Dirac cone requires an OP magnetic 
field B± or magnetization M±. On the other hand, when 
the SS has a finite OP spin component Sz, application of 
an IP magnetic field _B|| would also lead to a gap opening 
due to the coupling between By and Sz [!]• This way 
of gapping would be useful for achieving a massive Dirac 
state since it can essentially avoid the complicated orbital 
effects which generally couple to the IP k's and spins . 
Hence our observation of the OP spin component assures 
a promising direction for realizing novel quantum state 
and device applications that require a gapped SS. 



In conclusion, our SR-ARPES measurements of Bi2Te3 
provide direct and compelling evidence for the OP spin 
component on the gapless SS, by simultaneously deter- 
mining the IP and OP spin textures at various k points in 
the Brillouin zone. Moreover, we found that a stronger 
hexagonal deformation of the FS leads to a larger OP 
spin component, by comparing the observed FS and the 
SR-ARPES data with those of TlBiSe2. Intriguingly, the 
OP polarization is found to be much smaller than that ex- 
pected from existing theory in both Bi2Te3 and TlBiSe2, 
which poses an interesting question for future studies. 

We thank T. Arakane for the development of an in- 
situ sample rotation system for the SR-ARPES measure- 
ments and Zhi Ren for the preliminary crystal growth of 
Bi2Te3. This work was supported by JSPS (KAKENHI 
19674002 and the Next-Generation World-Leading Re- 
searchers Program), JST-CREST, MEXT of Japan (In- 
novative Area "Topological Quantum Phenomena"), and 
AFOSR (AOARD 10-4103). 

Note added. Recently, we became aware of a related 
SR-ARPES study on more highly electron-doped Bi2Te3 
22| . The reported Pz value (0.3) is larger than that of 
the present study (0.15). This difference is consistent 
with the k-p theory and is attributed to the difference 
in the magnitude of the warping factor w [2^ between 
the two cases (1.1 in the present sample, as opposed to 
1.6 in Ref. [H). 
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